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, Orbit functions of a simple Lie group/Lie algebra L consist of exponential functions summed 

' up over the Weyl group of L. They are labeled by the highest weights of irreducible finite 

O ^. dimensional representations of L. They are of three types: C-, S- and £'-functions. Orbit 

D ■ functions of the Lie algebras A„, or equivalently, of the Lie group SU(n+ 1), are considered. 

\ First, orbit functions in two different bases - one orthonormal, the other given by the simple 

' roots of SU(n) - are written using the isomorphism of the permutation group of n elements 

CN \ and the Weyl group of SU(ri). 

Secondly, it is demonstrated that there is a one-to-one correspondence between classical 
Chebyshev polynomials of the first and second kind, and C- and S'-functions of the simple 
Lie group SU(2). 

It is then shown that the well-known orbit functions of SU(n) are straightforward gener- 
■ alizations of Chebyshev polynomials to n — 1 variables. Properties of the orbit functions 

\ provide a wealth of properties of the polynomials. 

Finally, multivariate exponential functions are considered, and their connection with orbit 
functions of SU(n) is established. 

(N ' 
> 

CNj ! 1 Introduction 
0\ _ 

, The history of the Chebyshev polynomials dates back over a century. Their properties and 

\ applications have been considered in many papers. We refer to [19, 20] as a basic reference. 

' Studies of polynomials in more than one variable were undertaken by several authors, namely 

. [2-4,13,15,21,22]. Of these, none follow the path we have laid down here. 

In this paper, we demonstrate that the classical Chebyshev polynomials in one variable are 
naturally associated with the action of the Weyl group of SU(2), or equivalently with the action 
of the Weyl group IF(Ai) of the simple Lie algebra of type Ai. The association is so simple that 
5^ I it has been ignored so far. However, by making the cornerstone of our rederivation of 

Chebyshev polynomials, we have gained insight into the structure of the theory of polynomials. 
In particular, the generalization of Chebyshev polynomials to any number of variables was a 
straightforward task. It is based on the Weyl group W{An), where n < oo. This only recently 
became possible, after the orbit functions of simple Lie algebras were introduced as useful special 
functions [18] and studied in great detail and generality [8,9,11]. 

We proceed in three steps. In Section 2, we exploit the isomorphism of the group of per- 
mutations of n -|- 1 elements S and the Weyl group of SU(n + 1), or equivalently of An, and 
define the orbit functions of A„. This opens the possibility to write the orbit functions in two 
rather different bases, the orthnormal basis, and the basis determined by the simple roots of 
An, which considerably alters the appearance of the orbit functions. In the paper, we use the 
non-orthogonal basis because of its direct generalization to simple Lie algebras of other types 
than An- 

In Section 3 we consider classical Chebyshev polynomials of the first and second kind, and 
compare them with the C- and S'-orbit functions of Ai. We show that polynomials of the 
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first kind are in one-to-one correspondence with C-functions. Polynomials of the second kind 
coincide with the appropriate S'-function divided by the unique lowest non-trivial S'-function. 
We point out that polynomials of the second kind can be identified as irreducible characters 
of finite dimensional representations of SU(2). Useful properties of Chebyshev polynomials can 
undoubtedly be traced to that identification, because the fundamental object of representation 
theory of semisimple Lie groups/algebras is character. In principle, all one needs to know about 
an irreducible finite dimensional representation can be deduced from its character. An important 
aspect of this conclusion is that characters are known and uniformly described for all simple Lie 
groups / algebras. 

In Section 4 we provide details of the recursive procedure from which the analog of the 
trigonometric form of Chebyshev polynomials in n variables can be found. Thus there are n 
generic recursion relations for An, having at least n + 2 terms, and at most ^[(„+iJ/2]) + ^ 
terms. Irreducible polynomials are divided into n+1 exclusive classes with the property that 
monomials within one irreducible polynomial belong to the same class. This follows directly 
from the recognition of the presence and properties of the underlying Lie algebra. 

In subsection 4.2, the simple substitution z = e^'^*^, x G W^, is used in orbit functions 
to form analogs of Chebyshev polynomials in n variables in their non-trigonometric form. It is 
shown that, in the case of 2 variables, our polynomials coincide with those of Koorwider [13] (III), 
although the approach and terminology could not be more different, ours being purely algebraic, 
having originated in Lie theory. 

In Section 5, we present the orbit functions of An disguised as polynomials built from multi- 
variate orbit functions of the symmetric group. In Section 2, such a possibility is described in 
terms of related bases, one orthonormal (symmetric group), the other non-orthogonal (simple 
roots of An and their dual w-basis). Both forms of the same polynomials appear rather different 
but may prove useful in different situations. 

The last section contains a few comments and some questions related to the subject of this 
paper that we find intriguing. 

2 Preliminaries 

This section is intended to fix notation and terminology. We also briefly recall some facts about 
S„+i and An, dwelling particularly on various bases in and M". In Section 2.3, we identify 

elementary reflections that generate the An Weyl group W, with the permutation of two adjacent 
objects in an ordered set of n -I- 1 objects. And, finally, we present some standard definitions 
and properties of orbit functions. 

2.1 Permutation group S^+i 

The group S„+i of order {n + 1)! transforms the ordered number set [li,l2, ■ ■ ■ ,ln,ln+i] by 
permuting the numbers. 

We introduce an orthonormal basis in the real Euclidean space M"^"*"^, 

eieW+\ {ei,ej) = dij, l<i,j<n + l, (1) 

and use the l^s as the coordinates of a point fi in the e-basis: 

n+1 

k=i 

The group S^+i permutes the coordinates 1^ of /i, thus generating other points from it. The 
set of all distinct points, obtained by application of S„+i to fi, is called the orbit of S„+i. We 
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denote an orbit by W\, where A is a unique point of the orbit, such that 

h > h > ■ ■ ■ > In ^ In+l ■ 

If there is no pair of equal l^'s in A, the orbit Wx consists of (n + 1)! points. 

Further on, we wih only consider points f^ from the n-dimensional subspace TL C M""*"^ defined 
by the equation 

n+1 

J2lk = 0. (2) 

fe=i 



2.2 Lie algebra A„ 

Let us recall basic properties of the simple Lie algebra An of the compact Lie group SU(n + 1). 
Consider the general value (1 < n < oo) of the rank. The Coxeter-Dynkin diagram, Cartan 
matrix £, and inverse Cartan matrix of An are as follows: 



ooo- 

a.1 a.2 03 



O-n—l CX-n 



2-1000 
-1 2-1 
0-1 2-1 



0\ 
\ 




... -1 2 -1 / 
... -1 2/ 



n + 1 



l-n l-(n-l) l-(n-2) l-(n-3) 

l-(n-l) 2-(n-l) 2-(n-2) 2-(n-3) 

l-(n-2) 2-(n-2) 3-(n-2) 3-(n-3) 

l-(n-3) 2-(n-3) 3-(n-3) 4-(n-3) 



V 



1-3 
1-2 
1-1 



2-3 
2-2 
2-1 



3-3 4-3 
3-2 4-2 
31 41 



1- 3 1-2 

2- 3 2-2 

3- 3 3-2 

4- 3 4-2 



1- 1 

2- 1 
31 
4-1 



(n-2)-3 {n-2)-2 {n-2)-l 
(n-2)-2 (n-l)-2 (n-l)-l 
(n-2)-l (n-l)-l n-1 / 



The simple roots ttj, 1 < i < n oi An form a basis (a-basis) of a real Euclidean space M". 
We choose them in H: 

Oii = Ci- ej+i, i = 1, . . . ,n. 

This choice fixes the lengths and relative angles of the simple roots. Their length is equal to a/2 
with relative angles between and a^+i <:k < n—\) equal to and ^ for any other pair. 

In addition to e- and ct-bases, we introduce the a;-basis as the Z-dual basis to the simple 
roots ai'. 

{ai,Uj)=5ij, l<i,j<n. 

It is also a basis in the subspace H C M"+^ (see (2)). The bases a and oj are related by the 
Cartan matrix: 

a = €10, CO = €-~^a. 



Throughout the paper, we use X E H- Here, we fix the notation for its coordinates relative 
to the e- and w-bases: 

n+1 n n+1 

A = ^ IjCj =: (Zi, . . . , ln+i)e = ^ AjWi =: (Ai, . . . , Xn)uj, ^ k = 0. 

j=l 1=1 i=l 
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Consider a point X G H with coordinates Ij and Aj in the e- and w-bases, respectively. Using 
a = €uj, i.e. coi = Ylk=ii^~^)ik(^kj we obtain the relations between Aj and Iji 



k=l k=l 

/,=Ai(£r--G:ri-i)+A2(C2-/-£2i-i)+- • ■+KK]-^nU), J = 2, . . . ,n. 

or explicitly, 

Aj = — li+i, i = 1,2, ... ,n. 
The inverse formulas are much more complicated 
l = AX, 

where I = {li, . . . , In+i), A = (Ai, . . . , A„), and A is the (n+1) x n matrix: 



(3) 
(4) 



n+l 



/ n ft— 1 n—2 ■■■ 
I -1 n-1 n-2 ••• 
-1 -2 n-2 ••• 



-1 -2 -3 - -(n-1) 1 
\-l -2 -3 ••• -(n-l)-n/ 



2.3 The Weyl group of A„ 

The Weyl group W{An) of order (n + 1)! acts in H by permuting coordinates in the e-basis, 
i.e. as the group S„+i. Indeed, let ri, 1 < i < n he the generating elements of W{An), 
i.e, reflections with respect to the hyperplanes perpendicular to ai and passing through the 

n+l 

origin. Let x = ^k^k = {xi,X2, . . . ,Xn+i)e and (•, •) denote the inner product. We then have 

k=i 

the reflection by rf. 



riX = X 



-{x,ai)ai = {xi,X2, ■ ■ .,Xn+i)e - {xi - Xj+i)(ej - ej+i) 



(xi , . . . , Xi—\, Xi^\ , Xi, Xi-^-2 , • • • , 2^n+l )e- 



(5) 



Such transpositions generate the full permutation group S„+i. Thus, W{An) is isomorphic to 
Sn+i, and the points of the orbit W\{Sn+i) and Wx{An) coincide. 



2.4 Definitions of orbit functions 

The notion of an orbit function in n variables depends essentially on the underlying semisimplc 
Lie group G of rank n. In our case, G = SU(n+ 1) (equivalcntly, Lie algebra An). Let the basis 
of the simple roots be denoted by a, and the basis of fundamental weights by co. 
The weight lattice P is formed by all integer linear combinations of the w-basis, 

P = Zwi + ZCJ2 + • • • + '^^n- 

In the weight lattice P, we define the cone of dominant weights P~^ and its subset of strictly 
dominant weights P"'""'" 

P D P+ = Z- V + • • • + Z-°w„ D P++ = Z> V + • • • + Z>°a;„. 

Hereafter, C W denotes the even subgroup of the Weyl group formed by an even number 
of reflections that generate W . W\ and are the corresponding group orbits of a point A G M"". 

We also introduce the notion of fundamental region F(G) C M". For A^ the fundamental 
region F is the convex hull of the vertices {0, (jJ\,U2, ■ ■ ■ , Wn}- 
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Definition 1. The C orbit function C\{x), A G is defined as 



Cx{x):= e^'^^^^'^), xGM". (6) 

Definition 2. Tlie 6" orbit function S\{x), X G P"*""*" is defined as 

Sx{x):= (_i)P(M)e27ri(/.,x)^ xGR", (7) 

where p(^) is the number of reflections necessary to obtain fi from A. Of course the same can 
be obtained by different successions of reflections, but all routes from A to /i will have a length 
of the same parity, and thus the salient detail given by p{ijl), in the context of an 5- function, is 
meaningful and unchanging. 

Definition 3. We define E orbit function Ex{x), A G P*^ as 

Ex{x):= Y e^^'^i'^''\ xGM", (8) 
n&wliG) 

where := P'^ U r^P"^ and rj is a reflection from W. 

If we always suppose that X,iJ.eP are given in the w-basis, and a; G is given in the a 

n n n 

basis, namely A = ^ ^j^jj A* = X] t^j^j^ ^j^l^j ^ ^ x = ^ Xjaj, xj G M, then the orbit 
functions of have the following forms 

n 

Cx{x)=Y^ -1 = ^ JJe^-^^-^ (9) 
ii&Wx iJ.eWx j=i 

n 

Sx{x)= Y^-'^y^^^^ = Y{-^f^^^\[f^^''''''\ (10) 

n 

. 27ri Y H^3 « ^ J^. 

Ex{x)=J2^ = ^ JJe^-W. (11) 

new^ new^ j=i 

2.5 Some properties of orbit functions 

For S functions, the number of summands is always equal to the size of the Weyl group. Note 
that in the 1-dimensional case, C-, S- and ^'-functions are respectively a cosine, a sine and an 
exponential functions up to the constant. 

All three families of orbit functions are based on semisimple Lie algebras. The number of 
variables coincides with the rank of the Lie algebra. In general, C-, S- and E- functions are finite 
sums of exponential functions. Therefore they are continuous and have continuous derivatives 
of all orders in R". 

The S'-functions are antisymmetric with respect to the (n— l)-dimensional boundary of F. 
Hence they are zero on the boundary of F. The C-functions are symmetric with respect to the 
{n — l)-dimensional boundary of F. Their normal derivative at the boundary is equal to zero 
(because the normal derivative of a C-function is an 5*- function) . 

For simple Lie algebras of any type, the functions Cx{x), Ex{x) and Sx{x) are eigenfunctions 
of the appropriate Laplace operator. The Laplace operator has the same eigenvalues on every 
exponential function summand of an orbit function with eigenvalue — 47r(A, A). 
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2.5.1 Orthogonality 

For any two complex squared integrable functions (j){x) and ipix) defined on the fundamental 
region F, we define a continuous scalar product 

{4>ix),^l;ix)) := J </.(x)VK^dx. (12) 

F 

Here, integration is carried out with respect to the Euclidean measure, the bar means complex 
conjugation and x E F, where F is the fundamental region of either W or (note that the 
fundamental region of is F^ = F D riF, where € W). 

Any pair of orbit functions from the same family is orthogonal on the corresponding funda- 
mental region with respect to the scalar product (12), namely 

{C\ix),Cx'ix)) = \Wx\-\F\-Sxy, (13) 
{Sx{x),Sx'ix)) = \W\-\F\-Sxx', (14) 
{Ex{x),Ey{x)) = \W^\ ■ \F^\ ■ 6xx', (15) 

where 6xx' is the Kronecker delta, [l^/'l is the order of the Weyl group, \Wx\ and \W^\ are the 
sizes of the Weyl group orbits (the number of distinct points in the orbit), and \F\ and \F^\ are 
volumes of fundamental regions. The volume |F| was calculated in [6]. 

Proof. Proof of the relations (13,14,15) follows from the orthogonality of the usual exponential 
functions and from the fact that a given weight /i G P belongs to precisely one orbit function. □ 

The families of C-, 5*- and ^'-functions are complete on the fundamental domain. The com- 
pleteness of these systems follows from the completeness of the system of exponential functions; 
i.e., there does not exist a function (p{x), such that {(p{x) , (p{x)) > 0, and at the same time 
{(l){x),ip{x)) = for all functions iIj{x) from the same system. 



2.5.2 Orbit functions of acting in 



pn+l 



Relations (4) allow us to rewrite variables A and x in an orbit function in the e-basis. Therefore 
we can obtain the C-, S- and E- functions acting in ]R"+^ 

Cx{x)= e^'^^^^W,^)^ (16) 

Cx{x)= J2 (sgns)e2-^(^W'^), (17) 

SSSn+l 

Ex{x)= Yl e^^'^'^^^'^^^ (18) 
seAit„+i 

where (• , •) is a scalar product in M"+^, sgns is the permutation sign, and Alt„+i is the alter- 
nating group acting on an (n -|- l)-tuple of numbers. Note that variables x and A are in the 
hyperplane TC. 

Using the identity (A, rjx) = (r^ A, x) for the reflection r^, f = 1, . . . , n, it can be verified that 
Cxinx) = Crix{x) = Cx{x), and S^xix) = Sxinx) = -Sx{x). (19) 

Note that it is easy to see for generic points that Ex{x) = \ ^Ca(x) + Sx{x)^ , and from the 
relations (19), we obtain 

Er,x{x) = Exinx) = \ (Cx{x) - Sx{x)) = Ex{x). (20) 

A number of other properties of orbit functions are presented in [8,9,11]. 
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3 Orbit functions and Chebyshev polynomials 

We recall known properties of Chebyshev polynomials [19] in order to be subsequently able to 
make an unambiguous comparison between them and the appropriate orbit functions. 

3.1 Classical Chebyshev polynomials 

Chebyshev polynomials are orthogonal polynomials which are usually defined recursively. One 
distinguishes between Chebyshev polynomials of the first kind T„: 

ro(x) = l, Ti{x)=x, r„+i(x) = 2xT„-T„_i, (21) 
hence T2{x) = 2x^-1, Ts{x) = Ax^ - 3x, . . . (22) 

and Chebyshev polynomials of the second kind Un'- 

Uoix) = l, Uiix) = 2x, Un+i{x) = 2xUn - Un-l, (23) 
in particular C/2(x) =4x^-1, U^lx) = 8x^ - 4x, etc. (24) 

The polynomials T„ and Un are of degree n in the variable x. All terms in a polynomial have 
the parity of n. The coefficient of the leading term of r„ is 2**"^ and 2** for f/„, n = 1, 2, 3, . . . . 

The roots of the Chebyshev polynomials of the first kind are widely used as nodes for poly- 
nomial interpolation in approximation theory. The Chebyshev polynomials are a special case of 
Jacobi polynomials. They are orthogonal with the following weight functions: 



1 



' 0, 




m, 


< TT, 


n = 


m = 0, 


TT 

^ 2 ' 


n = 


m 7^ 0, 




n + 


m, 


I, 2' 


n = 


m. 



(25) 



Vl-xWn{x)Um(.x)dx=\ "Jj^ (26) 



There are other useful relations between Chebyshev polynomials of the first and second kind. 
^Tn{x) = nUn-i{x), n= 1,2,3,... (27) 

Tnix) = \iUn{x) - Un-2{x)), n = 2, 3, . . . (28) 

Tn+i{x) = xTn{x) - {1 - X^)Un-l, n = 1,2,3,... (29) 
Tn{x) = Un{x) - xUn-i, n=l,2, 3, ... (30) 

3.1.1 Trigonometric form of Chebyshev polynomials 

Using trigonometric variable x = cos y, polynomials of the first kind become 

T„(x) = T„(cosy) = cos(ny), n = 0,1,2,... (31) 
and polynomials of the second kind are written as 

Un{x) = Un{cosy) = ^'''^^''.^^^y\ n = 0,1,2,... (32) 

smy 

For example, the first few lowest polynomials are 

To{x) = To(cosy) = cos(Oy) = 1, Ti{x) = Ti{cosy) = cos(y) = x, 
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T2{x) = r2(cos y) = cos(2j/) = cos^ y — sin^ y = 2cos^y — 1 = — 1; 

Uo(x) = C/o(cosy) = = 1, Ui{x) = C/i(cosy) = . = 2cosy = 2x, 
sm y sm y 

, , , , sin(3u) sm(2y) cos y + sin'(/cos(2v) ,9 , . 9 , 

i72(x) = U2(cosy) = — ^ = — ^-^ = 4cos^y - 1 = 4x^ - 1. 

sm y sm y 

3.2 Orbit functions of and Chebyshev polynomials 

Let us consider the orbit functions of one variable. There is only one simple Lie algebra of 
rank 1, namely Ai. Our aim is to build the recursion relations in a way that generalizes to 
higher rank groups, unlike the standard relations of the classical theory presented above. 



3.2.1 Orbit functions of Ai and trigonometric form of T„ and Un 

The orbit of A = muJi has two points for m 7^ 0, namely W\ = {{m), (— m)}. The orbit of A = 
has just one point, Wq = {0}. 

One-dimensional orbit functions have the form (see (9), (10), (11)) 

C^(x) = e2^^'"^+e-2'^^'"^=2 cos(27rma:)=2 cos(my), where y = 2Trx, m e (33) 
5;,(a;) = e2~-e-2™=2zsin(27rmx)=2isin(my), for m e Z>°; (34) 
E^{x) = e^'^^'"^ = y"", where y = e^^^^^, m e Z. (35) 

Prom (33) and (31) it directly follows that polynomials generated from Cm functions of Ai are 
doubled Chebyshev polynomials of the first kind for m = 0, 1, 2, ... . 

Analogously, from (34) and (32), it follows that polynomials "g^^ are Chebyshev polynomials 
Um of the second kind for m = 0, 1, 2, . . . . 

The polynomials generated from Em functions of Ai, form a standard monomial sequence 
y"^, m = 0,1,2 ... , which is the basis for the vector space of polynomials. 

C- and S-orbit functions are orthogonal on the interval F = [0, 1] (see (13) and (14)) what 
implies the orthogonality of the corresponding polynomials. 

Comparing the properties of one-dimensional orbit functions with properties of Chebyshev 
polynomials, we conclude that there is a one-to-one correspondence between the Chebyshev 
polynomials and the orbit functions. 



3.2.2 Orbit functions of Ai and their polynomial form 

In this subsection, we start a derivation of the Ai polynomials in a way which emphasizes the role 

of the Lie algebra and, more importantly, in a way that directly generalizes to simple Lie algebras 
of any rank n and any type, resulting in polynomials of n variables and of a new type for each 
algebra. In the present case of Ai , this leads us to a different normalization of the polynomials 
and their trigonometric variables than is common for classical Chebyshev polynomials. No new 
polynomials emerge than those equivalent to Chebyshev polynomials of the first and second 
kind. Insight is nevertheless gained into the structure of the problem, which, to us, turned out 
to be of considerable importance. We are inclined to consider the Chebyshev polynomials, in 
the form derived here, as the canonical polynomials. 

The underlying Lie algebra Ai is often denoted sZ(2,C) or su{2). In fact, this case is so 
simple that the presence of the Lie algebras has never been acknowledged. 

The orbit functions of Ai are of two types (33) and (34); in particular, Co (a;) = 2, and 
So{x) = for all x. 
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The simplest substitution of variables to transform the orbit functions into polynomials is 
y = e^'^*^, monomials in such a polynomial are y"^ and y~^. Instead, we introduce new ('trigono- 
metric') variables X and Y as follows: 

X ■= Ci(x)=e2'^^^+e-2'^^^=2 cos(27rx), (36) 
Y ■= 5i(a;)=e2''^^-e-2'^^^=2zsin(27rar). (37) 

We can now start to construct polynomials recursively in the degrees of X and Y, by calculating 
the products of the appropriate orbit functions. Omitting the dependence on x from the symbols, 
we have 

X2 = C2 + 2 =^ C2 = X2 - 2, 

XC2 =C^ + X =^ C3 = X3 - 3X, (38) 

XCm = Cm+l + C'm-l =^ C'm+l = XCm — Cm-1, m > 3. 

Therefore, we obtain the following recursive polynomial form of the C-functions 

Co = 2, Ci = X, C2 = X^-2, C3 = X^-3X, C4 = X'^-AX'^ + 2,.... (39) 

After the substitution z = ^X we have 

Co=2-l, Ci=2z, ^2=2(2^2-1), C3=2{4z^-3z), 04=2(8^^-8x^ + 1),.... 

Hence we conclude that Cm = 2Tm, for m = 0, 1, . . . . 



Remark 1. 

In our opinion, the normalization of orbit functions is also more 'natural' for the Chebyshev 

'2 



polynomials. For example, the equality C| = C4 + 2 does not hold for T2 and r4. 



Remark 2. 

Each Cm also can be written as a polynomial of degree m in X,Y and 5^-1. It suffices to 
consider the products YSm, e.g., C2 = Y^ + 2, C3 = YS2 + X, etc. Equating the polynomials 
obtained in such a way with the corresponding polynomials from (38), we obtain a trigonometric 
identity for each m. For example, wc find two ways to write C2, one from the product X^ and 
one from Y^. Equating the two, we get 

X^-Y^ = 4 ^ sin2(27rx)+cos^(27ra;) = 1 

because Y is defined in (37) to be purely imaginary. 

Just as the polynomials representing Cm were obtained above, it is possible to to find poly- 
nomial expressions for Sm for all m. 

Fundamental relations between the S- and C- orbit functions follow from the properties of 
the character Xmix) of the irreducible representation of Ai of dimension m + 1. 

The character can be written in two ways: as in the Weyl character formula and also as the 
sum of appropriate C-functions. Explicitly, we have the Ai character: 

Sm+i{x) ^ / X , ^ / X , , iC2ix) + 1 for m even, 

^i{X) IC3(x) + Ci(x) for m odd. 

Let us write down a few characters 

- SIM - 1' - SIR - <-l - ^, X2 - - 62 + Co - X - 1, 

X3 = m = C3 + Ci = - 2X, X4 = ||) = C4 + C2 + Co = - 3X2 + 1, . . . 
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Again, the substitution z = transforms these polynomials into the Chebyshev polynomials 
of the second kind = Um, m = 0, 1, . . . , indeed 

Remark 3. 

Note that in the character formula we used Cq = 1, while above (see (11) and (39)) we 
used Co = 2. It is just a question of normalization of orbit functions. For some applica- 
tions/calculations it is convenient to scale orbit functions of non-generic points on the factor 
equal to the order of the stabilizer of that point in the Weyl group W{Ai). 

4 Orbit functions of An and their polynomials 

This section proposes two approaches to constructing orthogonal polynomials of n variables 
based on orbit functions. The first comes from the decomposition of Weyl orbit products into 
sums of orbits. Its result is the analog of the trigonometric form of the Chebyshev polynomials. 
The second approach is the exponential substitution in [8] . 

4.1 Recursive construction 

Since the C- and S- functions are defined for of any rank ri = 1, 2, 3, it is natural to 
take C-functions and the ratio of iS-functions as multidimensional generalizations of Chebyshev 
polynomials of the first and second kinds respectively 

Tx{x) := Cxix), X e R", 

^a(x):=%3^, p = C^i+L^2+...+C^n = (l,l,...,l)c., XGM", 

where A is one of the dominant weights of Aji. 

The functions Tx and Ux can be constructed as polynomials using the recursive scheme 
proposed in Section 3.2.2. In the n-dimensional case of orbit functions of An, we start from the 
n orbit functions labeled by the fundamental weights, 

Xi:=C^,{x), X2:=C^,{x), X„ := C^Jx) , x e . 

By multiplying them and decomposing the products into the sum of orbit functions, we build 

the polynomials for any C- and S'-function. 

The generic recursion relations are found as the decomposition of the products Xi^.C(^ai,a2....,an) 
with 'sufficiently large' ai, 02, . . . , a„. Such a recursion relation has ( "^'^ ) +1 terms, where ( "j"^ ) 
is the size of the orbit of u>j . 

An efficient way to find the decompositions is to work with products of Weyl group orbits, 
rather than with orbit functions. Their decomposition has been studied, and many examples 
have been described in [5]. It is useful to be aware of the congruence class of each product, 
because all of the orbits in its decomposition necessarily belong to that class. The congruence 
number ^ of an orbit A of An, which is also the congruence number of the orbit functions Cx 
and Sx, specifies the class. It is calculated as follows, 

n 

01,02,. ..,a„)(^)) — #('S'(oi,02,...,o„)(3^)) — ^^kak mod (n + 1). (40) 

k=l 

In particular, each Xj, where j = 1, 2, . . . , n, is in its own congruence class. During the multi- 
plication, congruence numbers add up mod n + 1. 

Polynomials in two and three variables originating from orbit functions of the simple Lie 
algebras A2, C-z, G2, A3, Bs, and C3 are obtained in the forthcoming paper [17]. 



10 



4.2 Exponential substitution 

There is another approach to multivariate orthogonal polynomials, which is also based on orbit 
functions. Such polynomials can be constructed by the continuous and invertible change of 
variables 

yj = e^'''^i, x,gM, j = l,2,...,ri. (41) 

Consider an An orbit function Cx{x), S\{x) or E\{x), when A is given in the oj-basis and x is 

n 

given in the a-basis. Each of these functions consists of summands W ^^miijXj ^ where jij G Z 
are coordinates of an orbit point fi. Then the summand is transformed by (41) into a monomial 

n 

of the form Y\. Vj^ ■ It is convenient to label these polynomials by non-negative integer coordi- 

nates (mi, m2, . . . , m„) of the point A = miLUi + m2W2 + • • • rUntOn and to denote the polynomial 
obtained from the orbit function C\ as P^^_^ ^ ^ (analogously for S and E functions). Poly- 
nomials of two variables obtained from the orbit functions by the substitution (41) are already 
described in the literature [13], where they are derived from very different considerations. The 
detailed comparison is made in the following example. 

Example 1. Consider the A2 Weyl orbits of the lower weights (0, m)^^, (m, 0)^^ and the orbit 
of the generic point (mi, 7712)0), m,mi,m2 G Z^'^ 

^(o,m)(^2) = {(0,m), (-m,0), (m,-m)}, T4^(m,o)(^2) = {(m,0), (-m,m), (0, -m)}, 
W^(mi,m2)(^2) = {("T'i,m2)+, (-mi, mi+m2)", (mi+m2, -m2)", 

(-m2,-mi)~, (-mi-m2,mi)+, (m2, -mi-m2)^}. 

Suppose X = {xi, X2) is given in the a-basis, then the orbit functions assume the form 

C(o,o)(a;) = 1, C(o,„)(a;) = C^mfi){x) = e-2-mxi+g2.imxig-2™+g2™^ 

I „\ „27rimixi„27rjm2a;2 I „— 27rimixi„27ri(mi+m2)x2 I „27ri(mi+m2)xi„— 27rim2a;2 I 

*-"(mi,m2)V^/ — ^ ^ "T^ ^ -re ^ ' e -|- 

g— 27rjm2a;i ^—2mm\X2 _j_g— 27ri(mi +m2)a;i g27rimia;2 _^^mm2Xi ^—2ni{mi +m2)x2 (42) 

Q / \ 2mm\xi 2mm2X2 — 27rimia;i 27ri(mi+m2)a;2 2m{m\-\-m2)x]_ —2mm2X2 

^(mi,m2) V'''/ — e e e e e e 

g— 27rim23;ig— 27rimia;2_j_g— 27rj(mi+m2)a;i g27rimia;2_j_g27rim2a;ig— 27ri(mi+m2)a;2 

Using (41) we have the following corresponding polynomials 



^{0,0) — M),m — -fO,m — 2/1 +^1 2/2 +2/2 j 

^&,m2) = 2/l-2/r+2/r'"^2/r^^'"^^+2/!"^^"^^2/2 

yr"^^2/2-™H2/r^"^^+"^^^2/2™H2/r2/2"^"^^+"^^^ 
^(1,^2) = 2/r2/r-2/r^2/r^-''"^^-2/^^^"^^^2/2-'"^- 

2/r^2/2-"^H2/r^'"^+'"^^2/2"'^+2/r2/2"^'"^+'"'^- 

The polynomials e+ and e~ given in (2.6) of [13] (III) coincide with those in (42) whenever 
the correspondence a = 2ttxi, t = 21^X2 is set up. So, both the orbit functions polynomials of 
A2 and are orthogonal on the interior of Steiner's hypocycloid. 

It is noteworthy that the regular tessellation of the plane by equilateral triangles considered 
in [13] is the standard tiling of the weight lattice of A2. The fundamental region R of [13] 
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coincides with the fundamental region .^(^2) in our notations. The corresponding isometry 
group is the affinc Wcyl group of A2. 

Furthermore, continuing the comparison with the paper [13], we want to point out that orbit 
functions are eigenfunctions not only of the Laplace operator written in the appropriate basis, 
e.g. in a;-basis, the corresponding eigenvalues bring — 47r^(A, A), where A is the representative 
from the dominant Weyl chamber, which labels the orbit function. This property holds not only 
for the Lie algebra An and its Laplace operator, but also for the differential operators built from 
the elementary symmetric polynomials, see [8,9]. 

An independent approach to the polynomials in two variables is proposed in [22], and the 
generalization of classical Chebyshev polynomials to the case of several variables is also presented 
in [4]. A detailed comparison would be a major task because the results are not explicit and 
contain no examples of polynomials. 

5 Multivariate exponential functions 

In this section, we consider one more class of special functions, which, as it will be shown, 
are closely related to orbit functions of An- Such a relation allows us to view orbit functions 
in the orthonormal basis, and to represent them in the form of determinants and permanents. 
At the same time, we obtain the straightforward procedure for constructing polynomials from 
multivariate exponential functions. 

Definition 4. [12] For a fixed point A = {li,l2, ■ ■ ■ ,ln+i)e, such that h > h > ■ ■ ■ > In+i, 

n+l 

^ Zfc = 0, the symmetric multivariate exponential function oi x = {xi,X2, ■ ■ ■ ,Xn+i)e is 

k=l 

defined as follows 



:= det+ 



^2TTil2Xi ^2-iril2X2 ^2-Kil2X„+i 



^ g27riin+lXl ^1^iln+lX2 ^2mln+lXn+l j 



(43) 



Here, det"*" is calculated as a conventional determinant, except that all of its monomial terms 
are taken with positive sign. It is also called permanent [14] or antideterminant. 

It was shown in [12] that it suffices to consider D^{x) on the hyperplane a; G 7^ (see (2)). 
Furthermore, due to the following property of the permanent 

det+(ajj)^-=l = ^ ai,s(l)a2,s(2) • • • am,s(m) = XI '^s(l),lCts(2),2 • • • as{m),m 

we have 

_ ^ g27riZia;3(i) . . . g27riZ,„a;s(„+i) _ ^ ^2ni(X,s{x)) _ ^ ^2ni{s{\),x) _ 

Proposition 1. For all X,x G Ti. C M^^^, we have the following connection between the 
symmetric multivariate exponential functions in n + 1 variables, and C orbit functions of A^ 
D^(x) = kC\{x), where k = j^^^, \W\ and \Wx\ are sizes of the Weyl group and Weyl orbit 
respectively. In particular, for generic points, k = 1. 

Proof. Proof follows from the definitions of the functions C and (definitions 1 and 4 respec- 
tively) and properties of orbit functions formulated in Section 2.5.2. □ 
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Definition 5. [12] For a fixed point A = {hjhj - ■ ■ Jn+i)ej such that h > h > ■ ■ ■ > Wi) 

n+l 

^ Zfc = 0, the antisymmetric multivariate exponential function oi x = {xi,X2, ■ ■ ■ ,Xn+i)e G 
fe=i 

H is defined as follows 



D~{x) := det 



r,2Tvil2Xl 



r,2'IvihXn-\ 
^2-Kil2Xn 



+1 



V e 



2mln+lX\ p2mln+lX2 



,2l7iln+lXn+l 



(sgn s)e2"(^W'^), 

sGS„+i 



(44) 



where sgn is the permutation sign. 



Proposition 2. For all generic points X & H C M*^"*"^, we have the following connection 
D^{x) = Sx{x). 

The antisymmetric multivariate exponential functions D , and S orbit functions, equal zero 
for non-generic points. 

Proof. Proof directly follows from the definitions of functions S and D~ (definitions 2 and 5 
respectively), and properties of S functions formulated in Section 2.5.2. □ 

Definition 6. [7] The alternating multivariate exponential function D^'*(a;), for 
X = (xi, . . . , Xn+i)e, ^ = {h, ■ ■ ■ , ln+i)ej is defined as the function 



Df\x) := sdet 



^2iTilixi 
„2ml2X\ 



a2'KihX2 

J2ml2X2 



^2^r^ln+-iX^ „2niln+iX2 



rj2'JvillXn 
,2ml2Xn+l 



^27viln 



yiXn+1 



(45) 



where Alt^+i is the alternating group (even subgroup of S„+i) and 

, n+l 



27ri(A,ui(x)) 



sdet ^e^'^*'-''^*'^""'' := g27riiia::„(i)g27rii22;TO(2) . . . g27riZ„+ia;„(„+i) _ g 

W&Aitn+l ■UJSAltn+l 

Here, (A, x) denotes the scalar product in the (n + l)-dimensional Euclidean space. 

Note that Alt^ consists of even substitutions of Sm, and is usually denoted as Am', here we 
change the notation in order to avoid confusion with simple Lie algebra notations. 

It was shown in [7] that it is sufficient to consider the function L>^'*(x) on the hyper- 
plane H. xi + X2 + ■ ■ ■ + Xn+i = for A, such that h > h ^ h ^ ' " ^ J-n+i- 

Alternating multivariate exponential functions are obviously connected with symmetric and 
antisymmetric multivariate exponential functions. This connection is the same as that of the 
cosine and sine, with the exponential function of one variable D 



Alt/ 



HDt{x) + D-ix)). 



Proposition 3. For all generic points A G H C the following relation between the 

alternative multivariate exponential functions D^^^ and E-orbit functions of An holds true: 
D^'\x) = Ex{x). 

For non-generic points \, we have E\{x) = C\{x) and, therefore, E\(x) = kD^{x), where 
k- 1^1 



Proof. Proof directly follows from definitions 3 and 6, from the relation E = ^{C-\-S), and from 
the properties of orbit functions formulated in Section 2.5.2. □ 



13 



6 Concluding remarks 



1. Consequences of the identification of W-invariant orbit functions of compact simple Lie 
groups and multivariable Chebyshev polynomials merit further exploitation. It is conceiv- 
able that Lie theory may become a backbone of a segment of the theory of orthogonal 
polynomials of many variables. 

Some of the properties of orbit functions translate readily into properties of Chebyshev 
polynomials of many variables. However there are other properties whose discovery from 
the theory of polynomials is difficult to imagine. As an example, consider the decompo- 
sition of the Chebyshev polynomial of the second kind into the sum of Chebyshev poly- 
nomials of the first kind. In one variable, it is a familiar problem that can be solved by 
elementary means. For two and more variables, the problem turns out to be equivalent to 
a more general question about representations of simple Lie groups. In general the coef- 
ficients of that sum are the dominant weight multiplocities. Again, simple specific cases 
can be worked out, but a sophisticated algorithm is required to deal with it in general [16]. 
In order to provide a solution for such a problem, extensive tables have been prepared [1] 
(see also references therein). 

2. Our approach to the derivation of multidimensional orthogonal polynomials hinges on the 
knowledge of appropriate recursion relations. The basic mathematical property underlying 
the existence of the recursion relation is the complete decomposability of products of the 
orbit functions. Numerous examples of the decompositions of products of orbit functions, 
involving also other Lie groups than SU(n), were shown elsewhere [8,9]. An equivalent 
problem is the decomposition of products of Weyl group orbits [5]. 

3. Possibility to discrctizc the polynomials is a consequence of the known discretization of 
orbit functions. For orbit functions it is a simpler problem, in that it is carried out in the 
real Euclidean space R". In principle, it carries over to the polynomials. But variables of 
the polynomials happen to be on the maximal torus of the underlying Lie group. Only 
in the case of Ai, the variables are real (the imaginary unit multiplying the 5-functions 
can be normalized away). For An with n > 1 the functions are complex valued. Practical 
aspects of discretization deserve to be thoroughly investigated. 

4. For simplicity of formulation, we insisted throughout this paper that the underlying Lie 
group be simple. The extension to compact semisimple Lie group and their Lie alge- 
bras is straightforward. Thus, orbit functions are products of orbit functions of simple 
constituents, and different types of orbit functions can be mixed. 

5. Polynomials formed from £'-functions by the same substitution of variables should be 
equally interesting once n > 1. We know of no analogs of such polynomials in the standard 
theory of polynomials with more than one variable. Intuitively, they would be formed as 
'halves' of Chebyshev polynomials although their domain of orthogonality is twice as large 
as that of Chebyshev polynomials [11]. 

6. Orbit functions have many other properties [8,9,11] that can now be rewritten as properties 
of Chebyshev polynomials. Let us point out just that they are eigenfunctions of appropriate 
Laplace operators with known eigenvalues. 

7. Notions of multivariate trigonometric functions [10] lead us to the idea of new, yet to 
be defined classes of W-orbit functions based on trigonometric sine and cosine functions, 
hence also to new types of polynomials. 
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8. Analogs of orbit functions of Weyl groups can be introduced also for the finite Coxeter 
groups that are not Weyl groups of a simple Lie algebra. Many of the properties of orbit 
functions extend to these cases. Only their orthogonality, continuous or discrete, has not 
been shown so far. 

9. Our choice of the n dimensional subspace Ti. in W^~^^ by requirement (2), is not the only 
possibility. A reasonable alternative appears to be setting In+i = (orthogonal projection 
on M"). 
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